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Abstract 

In this paper, properties of a homogeneous Bose gas with a Feshbach resonance are studied 
in the dilute region at zero temperature. The stationary state contains condensations of atoms 
and molecules. The ratio of the molecule density to the atom density is 7rna 3 . There are two 
types of excitations, molecular excitations and atomic excitations. Atomic excitations are gapless, 
consistent with the traditional theory of a dilute Bose gas. The molecular excitation energy is 



finite in the long wavelength limit as observed in recent experiments on 85 Rb. In addition, the 

00 

■^J - ' decay process of the condensate is studied. The coefficient of the three-body recombination rate is 

in 

about 140 times larger than that of a Bose gas without a Feshbach resonance, in reasonably good 

o , 

£^ . agreement with the experiment on 23 Na. 
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I. INTRODUCTION 

Feshbach resonance!!! has been observed in 23 Na and a few other alkali sys- 

tems. In these systems, two atoms can collide and fall into a molecular state with a different 
spin. The scattering length can be tuned by changing the energy difference between the 
atomic and molecular states through a magnetic field. In the 85 Rb system, the scatter- 
ing length was tuned to the negative region to study the critical particle number of the 
condensate 0| . In the positive region, an unusual processes of particle loss were observed ^]. 
Recently oscillations between atomic and molecular states were discovered 

One of the motivations to study Feshbach resonance is to create a Bose system with a 
large and positive scattering length a. In the traditional theory of a dilute Bose gas 



every physical quantity is expressed as a power series in the gas parameter VSnna 3 , with 
higher order terms coming from high order perturbations. But the perturbation breaks down 
when the density n is of the order of 1/a 3 . It is a challenge to find out what happens to the 
system beyond the dilute region. In recent experiments p], na 3 can be tuned as high as 0.5. 

However, even in the dilute region na 3 <C 1, the system with a Feshbach resonance is 
different from the system without a Feshbach resonance because of the molecular states. The 
observation of oscillations between the atomic and molecular states p indicates macroscopic 
occupation of molecular states. In fact, in addition to photoassociation III, Feshbach 
resonance is considered as another method to study the mixture of atomic and molecular 
condensates. 

In Ref. 12j, a two-channel model was constructed to study Feshbach resonance, with the 



Hamiltonian given by 

h 2 h 2 

H = — V^ f • + V0 f • V0 + V - 7(0V 2 + h.c.) + htffa (1) 

2m Am 2 

where is the molecular field operator and ip is the atomic field operator. The magnetic 
detuning energy is given by h. The coupling constant of atomic fields is given by go- The 
coupling constant between atomic and molecular states is given by 7. The sign of 7 is chosen 
to be positive by gauge transformation. 
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121 . 1131 . 114 I15| of this model have been focused on condensate proper- 



Most studies 



ties. When multiple modes were included, the condensates were found to be depleted over 



;[16] 



time 161. In Ref. 



171 ]. a Hartree-Fock method was used to take into account the contribution 
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from atom pairs with opposite momentum, and the dynamics of the correlation between 
the condensates and atom pairs was studied numerically. However some properties of this 
model such as the complete excitation spectrum and the decay process are still unclear. 
In this paper, we study these problems in the dilute region where the scattering length a 
is much bigger than the scattering length far off the resonance but much smaller than 
the interparticle distance, florae 1/ 'yfn. A simple saddle-point method is developed to 
handle this case. 



II. CONDENSATE PROPERTIES 



In the two-channel system, the numbers of atoms and molecules are not constants. It is 
easier to work in the grand canonical ensemble with the grand thermodynamical potential 
given by 

F = V + 20 V), (2) 

where /i is the chemical potential. At zero temperature in the dilute region, both atoms 
and molecules condense. The dominant contribution to the grand potential comes from the 
condensation, given by 

= f l^ol 4 - 7(0^o + ex.) + h\<f> \ 2 - MW 2 + 2|0 O | 2 ), (3) 

where ipo and 0o are the expectation values of the atomic and molecular fields. The values 
of ipQ and 0o are determined from the saddle-point equations given by 

" - (g \H 2 - 270o -^o = 0, (4) 



(/i-2/i)0*- 7 ^ = O. (5) 



In the region with positive scattering lengths, h < and fi > 0, the nontrivial solution 
of the above equations is given by 

where g e ff = go — ^l 2 /{h — 2/i). This solution is a minimum of the grand potential in the 
subspace of ipo and a maximum in the subspace of 0o because molecules have lower energy 

n 

than atoms. It is dynamically stable in the condensate space [12]. 



Compared to the traditional theory, the relation between the chemical potential and the 
atomic condensate density given by Eq.© remains the same, with g e $ now as the effective 
coupling constant of the atoms. The scattering length is proportional to the effective coupling 
constant, 

a = mg eS /(4iTh 2 ) = a (l - — ), (8) 

h — 2/i 

where the scattering length in the absence of the resonance is given by a$ = mg /(47rh 2 ), 
and the constant A is given by A = 27 2 /q - The scattering length in the form of Eq.fjHJ) 
have been observed in the experiments |2|, la]- 

In the dilute region, the interparticle distance is much larger than any other length 
scale. Since the chemical potential fi is a function of the atomic condensate density, the 
absolute value of the detuning energy \h\ is much larger than the chemical potential fi. As 
a consequence, the molecular condensate density is much less than the atomic condensate 
density, 

„ 7 2 ^o 2 < JL <<1 (q) 

These condensate properties are consistent with the findings in Ref.|12|. 

As shown in the next section, due to the coupling between atoms and molecules, the 
total molecule density is bigger than the molecular condensate density by a factor of 1/Z as 
given in Eq. ([20|) . because the molecular eigenstates also contains atom pairs. However, in 
the dilute region, the total molecule density is still much smaller than the atomic condensate 
density. The ratio of the molecule density to the atom density n m /n a is approximately given 
by 

_ „ _ „ ma . (10) 

III. EXCITATIONS AND QUANTUM DEPLETION 

In the dilute region, the system contains two types of excitations, the molecular excitation 
and the atomic excitation. The bare energy of a molecule is given by the detuning energy 
h. The coupling to the atomic states strongly renormalizes the molecule energy, and the 
molecular eigenstate becomes a superposition of the bare molecular state and pairing states 
of atoms. This renormalization can be easily demonstrated in the vacuum case where there 
are few particles. 
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In the vacuum, at zero temperature, the propagator of the molecular field is given by 

G m (q,tt) = - r 1 — n ^ (11) 

U — -f — h + 2ji — S m (q, il) + id 

and the propagator of the atomic field is given by 

G.(k,w) = (12) 

cj — e fc + ^ + id 

where E m (q, O) is the self energy of the molecular propagator generated by the coupling 
between atoms and molecules, 

E m (q, n) = 2* 7 2 J^-J 70iG a (k, c)G Q (q - k, Q - u) (13) 
7 7 (27r) 3 fi-e fc -e k _ q + 2/i + i(5' ^ ^ 

and e k = % 2 k 2 /(2m). 

There is an unphysical ultra-violet divergence in the integral in Eq. (jl4j) . As shown in 
many examples in the effective field theory, the bare parameters are often divergent and the 
physical parameters are always finite. The divergence from the perturbation can be removed 
by introducing a divergent counter term in the Hamiltonian. (In the context of a dilute 
Bose gas without a Feshbach resonance, see Ref. ^J) Here we regularize the divergence by 
introducing a divergent part of the bare molecule energy — 7 2 / <i 3 fc/[2(27r) 3 efc]. After the 
regularization, the self energy is given by 

r d 3 k 1 1 

Mq, «) = 27 2 j — — — ^ + -) (15) 



,.2 



= 5(?^-*- a (16) 

The pole of the molecular propagator in Eq. fjlljl is now given by 

n, = ^ - 2p - 1 (J^ - Ah - - T ^4) 2 . (17) 
q 2 P 4 I V 4vr 2 ft 6 2irh 3 ^ V ; 

Far off the resonance, the molecule energy is simply given by h+e q /2 — 2/j,; near the resonance 
m 3 7 4 /(7r 2 ^, 6 ) 3> 4|/i|, the molecule energy is given by 

e q n 4n 2 h 6 h 2 
n q ^^--2fi — . 18 

2 m 3 7 4 
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Near the resonance and close to the pole, the propagator can be approximated as 

G ^ = lT^i6' < 19) 

where the renormalization factor Z is given by 

" 4rf ^ (20) 



rwy 2 V m 



In this region, the scattering length is approximately given by 



a m ' , , 21 
and the molecular binding energy is approximately given by 

fi o « n- - ^ 22 

The relation between the binding energy of molecules and the scattering length is similar to 
that in the single-channel problem, 

h 2 

fio « 2 - ^ ( 23 ) 

except that there is a mean-field shift by 1\x. 

So far the results presented in this section Eq. (|18H23)) are derived for the vacuum case, 
valid near the resonance when the scattering length a is much larger than the scattering 
length off the resonance do- For a system with many particles, these results are still valid 
as long as the system is in the dilute region ao < a <C 1/ ^/n, because in this region the 
many-body contribution to the molecular energy from the condensate is proportional to the 
particle density n, much smaller than \Qq\. However, beyond the dilute region when na 3 > 1, 
the many-body effect may be dominant and the results for the dilute case may no longer 
be accurate. In a recent experiment [3], a shift to the molecular energy larger than 2/i was 
observed, which may be due to the many-body effect. 

If the molecule density is initially slightly away from the stationary value due to a per- 
turbation of the molecule energy in the action a<f>'<f>, the molecule density will oscillate. The 
oscillation can be calculated in the linear response theory. In the long-time t — » oo limit, it 
is given by 

< f (t)0(t) > -<P 2 « -ia<f>l[< <M^o(0) > -0o] (24) 
w -a$ |^G m (0,w)e-**/» (25) 

« 2i ^02 e -<not/» (26) 
\n\ 



The oscillation frequency is simply given by the molecular binding energy Qq because of 
the molecular condensation. It is similar to the single-channel problem where the phonon 
frequency coincides with the energy of the single particle excitation. In the experiment 
the oscillation frequency fits perfectly with the molecule energy calculated from a microscopic 
potential [20]. This result was also obtained in the context of photoassociation(2lj], in the 



Hartree-Fock approach 



17| . and in Ref. 



22]. 

In the dilute region, most of the particles are in the atomic condensate. The long- 
wavelength and low-energy excitations are atomic excitations, described by the gaussian 
fluctuation 

8F 2 = £[(e fc - + 2g^l - f^)4Vfc + - 70o)(V^lk + h.c.)\. (27) 

The molecular part of the fluctuation is not included because molecular excitations have a 
gap VLq as shown above. However molecular excitations contribute a diagonal self-energy to 
the atomic propagator through the coupling between atoms and molecules 

E a (k, w) = 2 7 V 2 G m (A;, w) » (28) 

which is already included in the r.h.s. of Eq. (J27|) . 

After diagonalizing Eq. (j27|) . we obtain a linear dispersed mode 



5F 2 = C + J2E k clc k , (29) 



k 

where Ck = Ukipk + Vktp_ h , u\ = |(^^ + 1) and v\ = u\ — 1. The excitation energy is given 
by 



E k = ^(6^ + 2/1), (30) 

and the energy of quantum depletion C is given by 

C 8 ,m\x s 



M- (31) 



**=/^rf = ^® 1 ( 32 ) 



The density of the quantum depletion can be easily calculated. 

d 3 k 2 1 ,m/i^3 

These results are consistent with the first order results in the traditional theory. To this 
order, /1 = 4Trh 2 a/m, any physical quantity can be expressed in terms of the density of the 
atomic condensate ipQ and the scattering length a . 



IV. THE THREE-BODY RECOMBINATION 



h J (2vr) 3 
2tt r d 3 k 



In the dilute region, the condensate state is only metastable because occupation of molec- 
ular states can lower the system energy Most of the decay is due to the three-body recom- 
bination. In the two-channel system, the picture of this process is rather simple. One pair 
of atoms with opposite momentum in the quantum depletion break up, and one of the them 
interacts with an atom in the condensate to form a molecule leaving the other atom of the 
original pair as an excited atom. The energy and momentum can be conserved in this pro- 
cess because molecules have negative energies. Mathematically the decay comes from the 
perturbation of the atom-molecule coupling term — 27^o0kV , k i n the Hamiltonian. 

The probability density of this three-body recombination process can be calculated by 
Fermi's golden rule 

2tt r d 3 k 

" ' ' | <k|2 7 ^ o 0j> k |G> \ 2 5(E k + Q k ) (33) 

h J (27r) 3 4 ^o^^(^ + ^), (34) 

where \G > is the ground state of 5F 2 , |k > is the ground state plus a molecule with 
wavevector k and an excited atom with wavevector — k. As defined earlier, v k is the coherence 
factor, and Z is the probability of finding a molecule in a molecular eigenstate given in 
Eq. (f2T)|) . In the dilute region near resonance, the decay rate is approximately given by 

^ 2ix f d 3 k (9 u 2 /77~~,m. 3 ,,3efc h 2 , , , 

r * t/^ 16 ^^^ 3 ^'^ (35) 

= 32v / 3vr 2 -^ 6 a 4 . (36) 
m 

In Fig.Q, the particle-loss rate for the 23 Na system is plotted and compared with the 
experimental data 3]. In the region 2ao < a < 10ao (906G to 906. 9G for the magnetic 
field), the theoretical result is fairly close to the experimental data, although the theoretical 
derivation so far is for the homogeneous case and the experimental system is inhomogeneous. 
In the single-channel case 2sj], the rate of the three-body recombination is also proportional 
to n 3 a 4 , but the coefficient is about 3.9, about 140 times smaller than 32\/37r 2 , and cannot 
be fitted to the experimental result on 23 Na. The increase of the coefficient reflects the fact 
that in the two-channel system there is a virtual process allowing two atoms to collide and 
fall into a molecular state, so the three-body recombination process is greater. However, in 
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rubidium system, the experimental results are rather intricate 5}. A decay rate with a weak 
density dependence was observed, which seems to suggest that the decay is dominated by a 
process independent of the particle density. 

In addition, the collision between atoms and molecules can also provide a three-body 
decay but its decay rate is much smaller [24]. The decay of the condensate is consistent 
the numerical result in Ref. that when multiple modes are included the condensate 
is depleted over time. The energy scale of the decay rate is comparable to the chemical 
potential /i when ip^a? is of the order of 0.01, 

fiT 

Iv^TT^a 3 . (37) 



It indicates that the condensate may not be stable enough beyond the dilute region. 



V. CONCLUSION 



In conclusion, we have obtained properties of a dilute Bose gas with a Feshbach resonance 
at zero temperature. The system contains both atomic and molecular condensates but the 
atom density is much higher. The molecular excitation has a energy gap and the atomic 
excitations are gapless. The low-energy properties are consistent with the results of the 
traditional theory. However, the condensed state is only metastable. The coupling between 
atomic and molecular states gives rise to a three-body recombination process. The decay 
rate is proportional to ip^a 4 as i n the single- channel case, but the coefficient is much bigger 
as observed in the Na experiment . 

It is tempting to go beyond the dilute region in the current approach. However, there are 
two expansion parameters in higher orders terms neglected in this approach, ip aoy/a from 
the atomic interaction and ipom 2, y 2 y/a/ (h 4 \h\) from the atom-molecule coupling. The latter 
is of the order of V'oV^ which is much larger than unity near the resonance, and causes 
the expansion to break down beyond the dilute region. 

a shift in the molecular energy was observed, which cannot be explained by any theories so 
far. It seems to suggest that there is a non-trivial many-body effect beyond the dilute region 
and more theoretical effort are needed. 
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FIG. 1: The coefficients of the particle- loss rate for Na system. N/N < n 2 > is plotted versus 
the magnetic field as in Fig. 2 in Ref. 2]. The straight line is the theoretical result computed from 
Eq. (|36[) : the squares and triangles are the experimental data taken at ramp speeds 0.13G/ms and 
0.31G/ms of the magnetic field^]. The scattering length as a function of the magnetic field is 
computed from Eq.Q with the parameters given in Ref.^]. 
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